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EFFECTS OF A MAGNETIC FIELD ON THE CONDUCTION HEAT TRANSFER AT THE 


STAGNATION POINT OF A PARTIALLY IONIZED ARGON GAS 
by James F. Schmidt 
Levis Research Center 


SUMMARY 

In a region of high electron-conduction heat flux (very high entry speeds) 
a magnetic field oriented parallel to the surface of the body can essentially 
eliminate the electron conduction heat flux and thereby reduce the total con- 
duction heat transfer to the stagnation-point region. This conduction heat 
transfer is calculated for the following argon-gas stagnation-point flow con- 
ditions: pressures of 1.0, 0.1, 0.01, and 0.01 atmosphere at temperatures of 

12 000°, 10 600°, 9500°, and 11 500° K, respectively. For the stagnation-point 
flow condition of 0.01 atmosphere and 11 500° K the maximum heat-transfer reduc- 
tion (55 percent) is obtained for a magnetic -field strength of 12 000 gauss. 

The nonequilibrium boundary layer may have an even greater potential than the 
equilibrium boundary layer for a large heat-transfer reduction. 


INTRODUCTION 

Consideration of atomspheric entry at velocities of 45 000 feet per second 
and greater must include the effect of the ionization phenomena on the body 
heating. The presence of electrons in the gas surrounding a reentry body 
greatly increases the thermal conductivity of the gas and thereby the conduc- 
tion heat flux to the body, depending upon the degree of ionization. Although 
the electrons have high mobility, which would increase the diffusion process 
of ion-electron pairs to the wall (where the large ionization energy can be 
released), the present report is concerned only with the conduction heat flux 
and thus gives a conservative estimate of the reduction possible by means of 
the properly oriented magnetic field. 

Many analyses are available for determining the effect on the stagnation- 
point heat transfer produced by the interaction of a magnetic field with an 
electrically conducting fluid in a laminar boundary layer (e.g., refs. 1 to 3). 
These analyses conclude that only by the effect of the magnetic field on the 
inviscid flow (decreased velocity gradient) can the heat flux to a body be 
reduced. However, Chapman and Cowling (ref. 4) have shown that, in general, 
for a plasma, a magnetic field oriented perpendicular to the temperature gra- 
dient can reduce the heat flux to a surface independent of any magnetic -field 
effects on the inviscid flow. Reference 4 (p. 336) assumes that the gas is 



composed of identically charged particles. In this report it is assumed that 
the results of reference 4 can be applied to the electron contribution to the 
conduction heat flux for a partially ionized gas in the stagnation-point region. 
In order to include the magnetic-field effects on the electron heat flux 
easily, the thermal conductivity and, thus, the conduction heat flux are sepa- 
rated into two parts, the neutral-atom and the electron thermal conductivity. 
Fay's mixture rule (ref. 5) for a partially ionized argon gas is used to cal- 
culate the conduction thermal conductivity of the neutral atoms and electrons. 
Although a partially ionized argon gas mixture is composed of ions as well as 
neutral atoms and electrons, the ion contribution to the conduction thermal 
conductivity is negligible until the gas becomes almost fully ionized (ref. 5). 
Argon (a simple monatomic gas) is selected for this study because its thermo- 
dynamic and transport properties have been explored experimentally as well as 
theoretically, even at very high temperatures. 

A two-dimensional partially ionized flow around a circular cylinder with a 
constant, uniform magnetic field (parallel to the cylinder axis and thereby 
perpendicular to the stagnation-point heat flux) is analyzed at the stagnation 
point. This flow situation can easily be duplicated by shock-tube experiments. 
Since the problem of the inviscid flow around the cylinder, including Hall and 
ion slip effects, is presently not solved, the correct velocity gradient param- 
eter at the edge of the boundary layer is not known. Therefore, solutions of 
the laminar -boundary -layer equations were obtained for several assumed values 
of the velocity-gradient parameter. 

The numerical calculations were performed by Kenneth A. Pew and are de- 
scribed in appendix C . 


ANALYSIS 

Assumptions 

Continuum flow in thermodynamic equilibrium is assumed for the two- 
dimensional boundary- layer flow. Additional assumptions are 

(1) The induced magnetic field is negligible. 

(2) The electrically conducting gas in the boundary layer maintains charge 
neutrality. 

(3) The diffusion of the reacting argon gas species is neglected. (This 
could be a first-order effect at the highest temperatures considered.) 

(4) There is negligible thermal radiation- energy absorption in the shock 
layer or boundary layer (ref. 6). 

(5) Thermal diffusion is neglibible (ref. 7). 

(6) Local similarity is assumed for solution of the boundary- layer equa- 
tions . 
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Basic Equations 

The general equations governing the flow of an electrically neutral con- 
ducting gas in a magnetic field are 

Continuity: 

div(pV) =0 (l) 

Momentum : 

DV 3 / N 

p = F - VP + div ¥ (2) 

Energy: 

pghaiva= i+i-i- 7 * E* (3) 


where p and P are the density and pressure, respectively, at any point in 
— * — > 

the gas, V is the velocity vector, F is the Lorentz force vector, equal to 
J “X B, B is the magnetic -field vector, 1 is the shear stress tensor, h is 
the static enthalpy, and q is the heat-flux vector. (Symbols are defined in 
appendix A. ) The effective electric field vector E from reference 8 is 


— > v — » — * 

E = E + V X B + 


HT 


m n 


e-im- 


m c 


l m 2 


e 0 m- 




VP 


pp 


m l m 2 


pfpg ™ ^o^-‘ 


vn n 


2 m l 


(4) 


where m is the particle mass, e is the particle charge, n is the number 
density, Vnq is the ion concentration gradient, the subscript 2 indicates an 
electron, and the subscript 1 indicates an ion. The applied electric field 
vector E is assumed herein to be zero, and the velocity vector V is essen- 
tially zero for the stagnation flow region. The last two terms in equation (4) 
are diffusion effects due to the partial pressure gradient and the concentra- 
tion gradient which are neglected in this analysis. Hence, the effective elec- 

— >-X- 

tric field vector E reduces to zero. 

For zero ion slip, the current density vector J, is given by reference 8 
as 


J = 




2,2 


B X E + 


ffl T 


m 




B(B * E*) 


-* . 2 2 
1 + CD T m 


(5) 


where o is the scalar electrical conductivity of the gas and oyr m is the 
Hall coefficient based on the momentum collision time (see appendix B). Since 
the effective electric field vector is zero in this study, the current density 
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vector (from eg. (5)) is also zero. The momentum and energy equations (2) 
and (3) in steady flow then reduce to 

Momentum : 


p(v • W) = -VP + div ¥ 


( 6 ) 


Energy: 


PV * Vh = -div q + § + V * VP 


(7) 


Scalar multiplication of the momentum equation (6) with V gives 

— > 

PV ' (V • W) = -V * VP + V * div i 


( 8 ) 


Substituting for V • VP from equation (8) in equation (7) yields the 
energy equation as 

PV • Vh = -div q + V • div ? + I -PV • (V • W) (9) 

If the usual boundary- layer assumptions (ref. 9) are made, div I, 

V • div t and £ are of the form 


div ¥ 
V • div 



I 



( 10 ) 

( 11 ) 

( 12 ) 


where the velocity u is in the x-direction around the circular cylinder. 

These equations are valid if the viscous stress is assumed not to be a function 
of the magnetic field. 


From appendix B, equation (B9), div q is 

_ _a_ f\ Sh\ _a_ / a e \ ah 

1V q " Sy \c p by) “ 3y \ Cp d + ^)) a y 


(13) 


where Cp is the specific heat of the gas mixture at constant pressure, A^ is 
the conduction thermal conductivity of the neutral atoms, Ag is the conduction 
thermal conductivity of the electrons, and oyr is the Hall coefficient. From 
appendix B, equation (B14), cot is 


COT 


15.20 B 

T A / 2 nl 


S EA 


a 


, 1 . 331X10“^ 
+ ip2 


/ \l/2 

3.943 + log[ — — ) 

\ 2/ J 


(14) 
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where T is the temperature in °K, a is the degree of ionization, ng is the 
electron number density in electrons per cubic centimeter, B is the magnetic 
field in gauss, and is the electron - neutral-atom collision cross 

section and is assumed constant ( 3XlO _d ® cm?) (see appendix B). 

After equation (10) is substituted in equation (6), the momentum equation 
in the x-direction becomes 


Pu g +Pv 


ot. 



(15) 


where 


££e to e 

ax “ “ p e u e dx 


(16) 


With the substitution of equations (ll) to (13) in equation (9), the 
energy equation in terms of total enthalpy (H = h + v? / 2) in the stagnation 
region reduces to 


Pu 


sh , „ sh a f \ sh\ , a ( a e \ 3h 

s + pv 5? - s?) + 5? L (1 + 


(17) 


Finally, the following system of laminar-boundary-layer equations (for the 
stagnation region) is to be solved: 


Pu 


SH 


Pu gii + 





Su 


o 


du e + a Su\ 

+ a ( sh\ 

^ \c p (l + o> 2 t 2 ) dy) 


(18a) 

(18b) 


(18c) 


Transformed Equations 


The following transformation equations (a slight modification of the trans- 
formation used by Lees (ref. 10)) for two-dimensional flow are used in the 
boundary- layer equations ( 18 ) : 


T) = 


p e< u e > l/2 j J s_ 


(p v Pw) 1 / 2 ( 2 x ) 1 / 2 J q 


ay 


X = X 


(19a) 

(19b) 
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where the subscript e indicates conditions at the edge of the boundary layer 
and the subscript w indicates wall conditions. The derivatives transform as 
follows : 

8 . ^ 1/Z 8 

S _ t d br\ 

^ dx ^ Sx 

and the stream function ijr, is defined so that 
Let 

t = (P w P w ) 1 / 2 (u e ) 1 / 2 (2x) 1 / 2 f (l9g) 

Then 

f 1 ( T) ) = ~ (I9h) 

u e 

Also, let 

g(Tj) = ^ (I9i) 

Applying the transformation equations (19) to the momentum and energy equa- 
tions (18b) and (18c) and reducing give 

Momentum : 

(Cf")‘ + ff" + {3^ - (f’) 2 + = 0 (20a) 

Energy: 

+ ^S') + fg , (l + |) = 0 (20b) 

where Pr w is the Prandtl number at the wall and (3, the velocity gradient 
parameter , is 


(l9e) 
(l9f ) 


(19c) 
( 19d) 
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( 21 ) 


2x 

u e cbc 


The density ratio 


P e /p is given by 


= 1 + q _T_ 
P 1 + a e T e 


( 22 ) 


The density-viscosity ratio C referred to wall conditions (T w = 300° K) 


is 


C = 


^W _j£_ 

P W P W (1 + a)T P w 


_ PM- _ 


(23) 


where the viscosity p is obtained from reference 11. The neutral-atom con- 
duction heat -flux parameter Dj± is 


Da = £_ C P^ M = _ h. 

Pw G p ( 1 + a )^ c p 


(24) 


Similarly, the electron-conduction heat-flux parameter Dg is 


]>E = 


£_ c p,w 
p w c p \r 
1 + od 2 t 2 


(25) 


The neutral-atom thermal conductivity 


me neuurax-auom unermax c 
ductivity Ag (from ref. 5) are 


a a = 


A* 


1+1.44 T 


0.16 a 

1 - a 


(26) 


A e - 




1 + 0.78X10 -4 yt 1 ~ ~ 

\ a 


(27) 


where A^, the unionized neutral-atom thermal conductivity, is approximated by 

A^ = 5 .8X10“ 7 T 3 / 4 (28) 

and A^,, the electron thermal conductivity for a fully ionized gas, is given by 
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( 29 ) 


l 4. 4X10~ 15 T 5 / 2 

1.2 4X10 4 T^/ 2 ~ 

■ (na) 1 /^ J 

The thermal conductivity at the wall 7^ is calculated from equation (26) 
with a = 0 and T = T^. 



The following boundary conditions are applied to the transformed momentum 
and energy equations (20): 


f(w) = 

f'(v) = 

= 0 

f'(ri) -* 1 

as 

T] -> oo 

f " ( T| ) 0 

as 

T) “> oo 

g(w) = g w 


g(Tl) -» 1 

as 

T) oo 

; ' ( r) ) -* o 

as 

T) oo 




> 


j 


(30) 


Surface Heat Transfer 


The conduction heat-transfer rate to a surface is 


1 



(31) 


Substituting the transformation equations (19c) and (l9i) and equation 
( 2l) in to equation (31) and reducing yield the heat-transfer rate at the stag- 
nation point: 


<ls = 



H e 4 


(32) 


The quantity (l/p)(du e /dx) is determined by the inviscid flow at the edge 
of the boundary layer. This inviscid flow is considered herein to be inde- 
pendent of the magnetic field, so that the ratio of heat-transfer rates with 
and without the magnetic field becomes 


= (S i>B=0 
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RESULTS AM) DISCUSSION 


Flow Property Variation Across Boundary Layer 

The thermodynamic properties used in evaluating the flow coefficients in 
equations (20) were obtained from reference 11 for the following four 
stagnation-point flow conditions : 

The flow coefficients in equations (20) were calcu- 
lated from equations (22) to (29) as a function of tem- 
pera- ture for a given pressure but are plotted (figs. 1 
and 2) as a function the enthalpy ratio g, which is 
one of the two dependent variables for the boundary- 
layer equations. Equations (20) were solved with an 
analog differential analyzer which gives results to four 
significant figures. The calculation method is described 
in appendix C. Equations (20) were also programed on a 
digital 7094 computer because the last stagnation-point 
flow condition was out of the range of accurate analog 
computer solutions. 


The variation of the density-viscosity ratio C and 
the density ratio P e /P (from eqs. (22) and (23)) as functions of the enthalpy 
ratio g, is shown in figure 1. The density-viscosity ratio falls to a very 
small value at the edge of the boundary layer (g = 1.0). Since this ratio is 
used only in the momentum equation in the present analysis, it has a negligible 
effect on the heat transfer. 


Pressure at 
edge of 
boundary- 
layer, 

P e> 

atm 

Temperature 
at edge of 
boundary 
layer, 

T e> 

°K 

1.0 

12 000 

.1 

10 600 

.01 

9 500 

.01 

11 500 


The variation of the sum of the neutral-atom and electron heat-flux param- 
eters (from eqs. (24) and (25)) as functions of the enthalpy ratio is shown in 
figure 2. The presence of the magnetic field is seen to reduce the total con- 
duction heat flux through the boundary layer. The sum of the neutral-atom and 
electron conduction heat-flux parameters can also be expressed by the identity 


d a + 


I>E 


= %(l + 



The factor 1 + I>e/D a shows more clearly the contribution of the electrons to 
the total heat flux (see fig. 3). Since 


I>E 



1 + <d 2 t 2 


from equations (24) and (25), the effect of the magnetic field in reducing the 
electron portion of the heat flux is shown in figure 3. The degree of ioniza- 
tion a is also plotted in figure 3 to show that the electron heat flux is 
significant when the degree of ionization becomes appreciable. 


The variation of the Hall coefficient tax and the thermal-conductivity 
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ratio A^/Aa in the boundary layer is shown in figure 4. At very small 
degrees of ionization the electron heat flux is negligible compared with the 
neutral-atom heat flux, and hence the magnetic field effect on the total heat 
flux is negligible even though the Hall coefficient becomes very large. For 
large degrees of ionization the thermal-conductivity ratio A^/A^ increases 
and a greater potential effect of the magnetic field on heat transfer is possi- 
ble. For comparable magnetic fields figure 4 shows a very large increase in 
the Hall coefficient as the pressure is reduced. This indicates that a signif- 
icant viscous heat-flux reduction is possible for an equilibrium boundary layer 
with low pressures (P e = 0.1 to 0.01 atm). A much larger heat-flux reduction 
appears possible for the low-pressure nonequilibrium boundary layer because in 
addition to large Hall coefficients the degree of ionization may also be larger 
and remain significant right up to the wall surface (ref. 12). 


Solutions of Boundary-Layer Equations 

Although the flow property variations in the boundary layer can provide 
qualitative effects on the heat transfer, the boundary- layer solutions for the 
highly nonlinear system of equations are required to predict quantitative re- 
sults. Figure 5 presents these solutions in the form of the heat-transfer 

reduction ratio ( q a )_./( q G )_, . as a function of the applied magnetic field and 
5 ° Jd=U 

stagnation flow condition. The solutions were calculated for values of the 
velocity-gradient parameter 0 ranging from 0 to 2.0. The single curve for 
each stagnation flow condition in figure 5 indicates that the ratio of heat 
transfer with and without a magnetic field is independent of 0. However, the 
actual heat-transfer rate, with or without a magnetic field, depends signifi- 
cantly on the value of 0 (see fig. 6). 

Figure 5 shows that the effect of the magnetic field on the heat-transfer 
ratio is accentuated at lower pressures and higher temperatures. For the two 
stagnation-point flow conditions at the same pressure (P e = 0.01 atm) but dif- 
ferent stagnation-point temperatures (9500° and 11 500° K), the degree of 
ionization is, respectively, 0.1177 and 0.6281. Therefore, the large differ- 
ence in the heat-transfer reduction ratio between these two flow conditions is 
attributed to the effect of ionization. For each stagnation-point flow condi- 
tion, the heat-transfer ratio levels off as the magnetic field is increased 
(see fig. 5). An explanation of this phenomenon is that, once the electron 
heat flux has been reduced to a negligible amount, the heat transfer is en- 
tirely due to the neutral atoms which are not affected by the magnetic field. 


Lewis Research Center, 

Rational Aeronautics and Space Administration, 
Cleveland, Ohio, December 3, 1965. 
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APPENDIX A 


SYMBOLS 



G 

first, second, and third time division multiplier inputs 
magnetic field 

density-viscosity ratio (eq. (23)) 
mean speed of electrons 

specific heat of gas mixture at constant pressure 
neutral-atom conduction heat-flux parameter (eq. (24)) 
electron conduction heat -flux parameter (eq. (25)) 

total substantial derivative, • + V • V 

effective electric field 
particle charge 
Lorentz force 

integral of velocity distribution across boundary layer (eq. (l9h)) 
total- enthalpy ratio, H/H e 
total enthalpy 
static enthalpy 

unit vectors in x- and y-directions 
electric current density 
Boltzman's constant 
particle mass 

number density of gas mixture 
ion concentration gradient 


static pressure 
defined by eq. (C4) 
defined by eg,. (C3) 

Prandtl number, Cpi-i/A 
heat-flux or heat-transfer rate 

p[p e /p - (f ' ) 2 ] 
position vector 

electron - neutral-atom collision cross section 
electron-ion collision cross section 
static temperature 
time 

velocity parallel to cylinder surface 

velocity gradient at edge of boundary layer 

total flow velocity 

velocity normal to cylinder surface 

coordinate parallel to cylinder surface; x = x 

coordinate normal to cylinder surface 

coordinate parallel to cylinder axis 

degree of ionization 

velocity-gradient parameter 

constant, 1.865 

similarity variable defined by eq. (19a) 
neutral-atom conduction thermal conductivity 
electron conduction thermal conductivity 
conduction thermal conductivity at wall 
coefficient of viscosity 



HSU 


I viscous dissipation function 

p mass density 

a scalar electrical conductivity 

t electron energy collision time 

T m electron momentum collision time 
shear stress 
i r stream function given by eq,. (I9g) 

co electron cyclotron frequency 

Subscripts : 

B with magnetic field 

B=0 without magnetic field 

e edge of boundary layer 

s stagnation flow 

w wall 

1 ion 

2 electron 
Superscripts : 

average 

vector 

— > 

tensor 

VV" derivatives with respect to r\ 



APPENDIX B 


CONDUCTION HEAT FHJX AND HALL COEFFICIENT 


Conduction Heat Flux 

Chapman and Cowling (ref. 4) have shown that a magnetic field (perpendic- 
ular to the conduction heat flux) will reduce the conduction heat flux for a 
gas composed of identically charged particles (electrons). This reduced elec- 
tron heat flux (ref. 4) is 


q z = -A. 


E 


ST cot ^ ST 

a? ' b~ a? 

1 + C0 2 T 2 


(Bl) 


For two-dimensional flow with a constant magnetic field in the 
z -direction, the electron heat-flux equation (Bl) becomes 


A 


*2 = - 


E 


1 + co 2 t 2 


/-? ST -? ST -» ST St\ 

( 1 S + 0 5? - SJ + 5?) 


(B2) 


The divergence of this electron conduction heat flux (eq. (B2)) gives 


S / A E ST\ 5 I A E bT' 
div *2 - "Sx U + m 2 x 2 57 J " ! + cd 2 t 2 Sy, 


AgCOT 


ST\ + 5 


AgCDT 


’ dx \ i + co 2 t 2 \ i + cd 2 t 2 dx 

Expanding the last two terms of equation (B3) gives 


(B3 ) 


a eP t St\ 5 ( a e cot St\ 


A 


E 


Sx Ip + ^,2^.2 by I by Ip + cd 2 t 2 ^ x / 


n j. 2 2 

1 + CO T 


s 

^ 57 


A 


E 


1 + C0 2 t 2 by Sx 


ST S(cot) _ T ST S 


A 


E 


A e 


^ dX Vl + C0 2 T 2 ; 


1 + co 2 t 2 


COT 



^e St S(cot) + mT St S | 

i + co 2 t 2 Sx Sx \i 


A 


E 


V 1 + co 2 t 2 ^ 


(B4) 
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where 


5 (cot) _ 5 (cot) BT 
By “ BT By 

B (cot ) _ B(cot) BT 
Bx ~ BT Bx 

B / a e \ B / a e \ Bt 

5y \i + co 2 t 2 /" + o^T 2 / 5 ^ 


A / A E \ = B / a e \ bt 

Bx yi + co 2 t 2 / ^\1 + co 2 t 2 /^ 


» 


Substituting equation (B5) into equation (B4) and reducing give 


B 


A 


'E 


JT , B 

COT =7— I + 


A 


E 


" Sx \i + co 2 t 2 ^y/ V 1 + 

Therefore, equation (B3) reduces to 


“ T i = 0 


. . d 

dlV *2 = “ 


Ae bt\ _b_ ( Ae bt] 

,2 t 2 Bx) ■ By U + C0 2 T 2 ByJ 


ll + OTT 


(B6) 


(B7) 


However , the first term in equation (B7) is very small compared with the second 
term. Hence, the divergence of the electron heat flux becomes 



/ a e Bt\ 
\i + co 2 t 2 


(B8) 


For a partially ionized gas the atoms still conduct heat, but play a 
diminishing roll as the degree of ionization increases. Equation (338 ) expresses 
the conduction heat flux of the electrons. Also, the magnetic field can only 
affect the charged particles (electrons) and not the neutral atoms. Therefore, 
the conduction heat flux (in terms of static enthalpy) is separated into two 
parts, the neutral-atom and the electron heat flux. The divergence of this 
conduction heat flux is 


div q = 


B_ 

Aa Bh\ 

b r 

Bh"| 

By 

\c p ByJ 

^y Lc p d + u 2 t 2 ) 

‘ByJ 


(B9) 
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Hall Coefficient 


The Hall coefficient required to evaluate the electron heat flux can be 
expressed by the following equation: 


e 2 B T m 
" m 2 £ e 


(BIO) 


where eg is the electron charge, mg is the electron mass, x m is the elec- 
tron collision time based on the momentum transfer, and £ e is a constant 
(1.865) for a fully ionized gas (ref. 13) that reduces the electron collision 
time based on the momentum transfer to the electron collision time based on 
the energy transfer. For a fully ionized gas is used in the partially 

ionized gas in the boundary layer for the following reasons : 

(1) At very small degrees of ionization the electron heat flux is negli- 
gible compared with the neutral-atom heat flux, and hence the Hall effect on 
the heat flux is negligible (no matter how large the Hall coefficient becomes). 

(2) At larger degrees of ionization (outer region of the boundary layer), 
the electron - neutral-atom collisions are negligible compared with the 
electron-ion collisions. Therefore, the partially ionized gas behaves like a 
fully ionized gas. 


The electron collision time x m based on the momentum transfer can be 
approximated (from ref. 14) by 



where #,,, the mean speed of the electrons, is 



(Bll) 


(B12) 


Although the electron - neutral-atom collision cross section is not 

constant with temperature (see ref. 15), a constant value (3XlCT-^ cm 2 ) is used, 
for the same reasons as those for using the fully ionized gas value of the con- 
stant £ e . The electron-ion collision cross section is 


S EI 


1.551X10" 5 

ij>2 


3.943 + log 



(B13) 
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APPENDIX C 


CALCULATION METHOD 
by Kenneth A. Pew 

The solution of equations (20) was obtained with the use of an EAI 16-31R 
analog computer. The computations were recorded with an EAI variplotter and a 
Brush recorder. 

The first method attempted was the standard analog approach of solving for 
the highest order derivatives : 



The parameters p g /p, C, D^, and Bg were curved fitted as a function of g. 
The problem was then complicated by also requiring the evaluation of C', D^, 
and Dg. Various curve fit routines were used on an IBM 7094 digital computer 
to construct analytic expressions for C, D A , and Dg. The analytic functions 
fo un d were differentiated with respect to g and the derivatives set into the 
analog circuit on manual diode function generators. These derivatives were 
multiplied by g' to develop C', D^_, and Dg, and integrated to form C, D A , 
and Dg. This process generated so much error in evaluating C', D^, and Dg 
that the results would be of little value . 

This large error was avoided by using a transformation (see ref. 5) which 
eliminated the need for C 1 , D^, and Dg . Let 


= g'(D A + %) 

(C3) 

0» f = Cf" 

The transformed equations then become 

(C4) 

+ 1 ) 

« " " D A + D E 

(C5) 

f (l + ^ m _ ^1 


(C6) 
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This made it possible simply to set the functions p e /p, C, D^, and Dg into 
the circuit on manual diode function generators (fig. 7). The circuit did 
require time scaling in order to slow the solution down by a factor of 5. 


A trial solution was started by estimating the initial conditions 
and f^. The computer would then be operated and the results inspected. If 
the boundary conditions were not met, g^ and/or f^ would be altered. This 
procedure was repeated until the boundary conditions (eq. (21)) were satisfac- 
torily met (±1 percent). Each computer operation required about 10 seconds. 
After the effects of the two initial conditions upon the system had become 
familiar, a solution could be obtained in approximately 5 minutes. 

The greatest error involved in the solution was in the setting of the 
manual diode function generators. This error is estimated to be less than 
2 percent. 
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Density-viscosity ratio, or density ratio, 



of ionization, 0.1179. 



(b) Pressure, 0. 1 atmosphere; 
temperature, 10 600° K; degree 
of ionization, 0.1163. 



(c) Pressure 0.01 atmosphere; (d) Pressure, 0.01 atmosphere; 

temperature, 9500° K; degree temperature, 11 300° K; degree 

of ionization, 0.1177. of ionization, 0.6281. 


Figure 1. - Variation of density -viscosity ratio and density ratio with en- 
thalpy ratio in boundary layer for various stagnation conditions. 
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(a) Pressure, 1.0 atmosphere; 

temperature, 12 000° K; degree 
of ionization, 0. 1179. 
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(b) Pressure, 0. 1 atmosphere; 
temperature, 10 600° K; degree 
of ionization, 0. 1163. 



(c) Pressure, 0.01 atmosphere; 
temperature, 9500° K; degree 
of ionization, 0. 1177. 


(d) Pressure, 0. 01 atmosphere; 
temperature, 11 500° K; degree 
of ionization, 0.6281. 


Figure 2. - Variation of sum of neutral -atom - and electron -conduction 
heat-flux parameters with enthalpy ratio in boundary layer for various 
stagnation conditions. 
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Enthalpy ratio, g 


(c) Pressure, 0. 01 atmosphere; (d) Pressure, 0. 01 atmosphere; 
temperature, 9500° K; degree temperature, 11 500° K; degree 

of ionization, 0. 1177. of ionization, 0.6281. 


Figure 3. - Variation of heat-flux -para meter ratio and degree of ionization 
with enthalpy ratio in boundary layer for various stagnation conditions. 
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Thermal -conductivity ratio, X E /X A 



of ionization, 0. 1179. 



temperature, 10 600° K; degree 
of ionization, 0. 1163. 



(c) Pressure, 0. 01 atmosphere; (d) Pressure, 0. 01 atmosphere; 

temperature, 9300° K ; degree of temperature, 11 500° K ; degree 

ionization, 0.1177. of ionization, 0.6281. 


Figure 4. - Variation of Hall coefficient and thermal -conductivity ratio 
with enthalpy ratio and magnetic-field strength for various stagnation 
conditions. 
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Derivative of enthalpy ratio at wall, 



Figure 5. - Effect of magnetic field on heat-transfer ratio 
for equilibrium laminar boundary layer at stagnation 
point. 



of ionization, 0.1179. 
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(W Pressure, 0. 1 at- 
mosphere; tempera- 
ture, 10 300° K; degree 
of ionization, 0. 1163. 



0 12 0 12 
Velocity-gradient parameter, p 

(c) Pressure, 0.01 at- (d) Pressure, 0.01 at- 
mosphere; tempera- mosphere; temperature, 

ture, 9500° K; degree 11 500° K; degree of 
of ionization, 0. 1177. ionization, 0. 6281. 


Figure 6. - Effect of velocity-gradient parameter on 
heat transfer without magnetic field for various 
stagnation conditions. 
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